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We consider the Scharnhorst effect (anomalous photon propagation in the Casimir vacuum) at 
oblique incidence, calculating both photon speed and polarization states as functions of angle. 
The analysis is performed in the framework of nonlinear electrodynamics and we show that many 
features of the situation can be extracted solely on the basis of symmetry considerations. Although 
birefringence is common in nonlinear electrodynamics it is not universal; in particular we verify that 
the Casimir vacuum is not birefringent at any incidence angle. On the other hand, group velocity is 
typically not equal to phase velocity, though the distinction vanishes for special directions or if one is 
only working to second order in the fine structure constant. We obtain an "effective metric" that is 
subtly different from previous results. The disagreement is due to the way that "polarization sums" 
are implemented in the extant literature, and we demonstrate that a fully consistent polarization 
sum must be implemented via a bootstrap procedure using the effective metric one is attempting 
to define. Furthermore, in the case of birefringence, we show that the polarization sum technique is 
intrinsically an approximation. 



PACS numbers: 12.20.Ds, 41.20.Jb; juant-ph/0010055 



I. INTRODUCTION 

In 1990 Scharnhorst demonstrated that the propaga- 
tion of light in the Casimir vacuum [H, ||, || is character- 
ized by an anomalous speed [Q. In fact photons propa- 
gating perpendicular to the plates travel at a speed cj_ 
which slightly exceeds the usual speed of light c in the 
Minkowski vacuum [| [| M . The propagation of pho- 
tons parallel to the plates instead occurs at the usual 
speed c || = c. 

Unfortunately this anomalous propagation is far 
too small to be experimentally detectable; the rela- 
tive modifications to the speed of light are of order 
10~ 2 a! 2 /(m e a) 4 , where a is the fine structure constant, 
a is the distance between the plates, and m c is the elec- 
tron mass. [We work in units such that % = c = 1. 
Greek and Latin indices run from to 3 and from 1 to 
3, respectively, r]^ = diag (— 1, +1, +1, +1) denotes the 
Minkowski metric] It is nevertheless an important point 
of principle that a quantum vacuum which is polarized by 
an external constraint can behave as a dispersive medium 
with refractive index n(oj) which remains less than unity 
to arbitrarily high frequency. Similar effects have been 
discovered in the case of quantum vacuum polarization 
in gravitational fields J8|, [Io|, [ll], and several other 
situations ]T^ |. 

It is to be emphasized that the Scharnhorst effect, 
albeit small, is of fundamental theoretical importance. 
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Though the present calculations are carried out in the 
"soft photon" approximation (wavelengths much larger 
than an electron Compton wavelength) there is an argu- 
ment based on dispersion relations which strongly sug- 
gests the actual signal speed is also modified (indeed, 
increased) Q] — the physics here is very different from 
that of the resonance-induced "apparent" superluminal 
velocities currently of experimental interest [p"3[ |l4f| ; here 
we are dealing with quantum polarization induced "true" 
superluminal velocities, albeit well outside the realm of 
present day experimental technique. 

A common feature of the gravitational analogs of the 
Scharnhorst effect is the presence of birefringence: Pho- 
tons with different polarizations propagate at different 
speeds. Similarly in generic situations nonlinear electro- 
dynamics often leads to birefringence |L5| [l6|, [TtJ . Never- 
theless the occurrence of birefringence in nonlinear elec- 
trodynamics is not universal. Indeed, in the Casimir vac- 
uum between parallel conducting plates the possibility of 
birefringence is tightly constrained: There is a residual 
(2+1) Lorentz invariance which prevents birefringence for 
photons that propagate parallel to the plates, and simi- 
larly there is an 0(2) rotational invariance that prevents 
birefringence for photons that propagate perpendicular 
to the plates. It is only for photons that propagate 
obliquely to the plates that there is even any possibility of 
birefringence. We shall be interested in a complete anal- 
ysis of both propagation speed and polarization states as 
a function of angle, and as a side effect report the (melan- 
choly) conclusion that birefringence is completely absent 
in the Casimir vacuum. As a consequence, an effective 
metric description is sufficient for completely specifying 
photon propagation in the Casimir vacuum — though it 
will in general not be sufficient for other more general 
vacuum states. The effective metric we find is subtly at 
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variance with previous results obtained by performing a 
"polarization sum" |ll| |l9) . The disagreement is due to 
the way that "polarization sums" are implemented in the 
extant literature, where they are defined in terms of the 
fiat Minkowski metric. We show that in the absence of 
birefringence a fully consistent polarization sum must be 
implemented in terms of a bootstrap procedure; using 
as input the effective metric one is attempting to define. 
Furthermore, in the presence of birefringence, we demon- 
strate that the polarization sum technique is intrinsically 
an approximation. Fortunately the differences first show 
up at order a 4 , and in almost all circumstances are com- 
pletely negligible. 

Throughout, we emphasize the use of symmetry argu- 
ments as a way of extracting general information that 
is (as much as possible) independent of the particular 
choice of Lagrangian. Finally, we have a few words to 
say concerning the utility of effective metric approaches 
in other contexts well beyond the Casimir vacuum. 



II. GENERAL FORMALISM 

Anomalous photon propagation can most easily be in- 
terpreted in terms of nonlinear electrodynamics. After 
integrating out virtual electron loops, the Maxwell La- 
grangian should be replaced (in the absence of bound- 
aries, at the one-loop level, and provided the distance 
scales defined by the field gradients and photon wave- 
length are much larger than the electron Compton wave- 
length) by Schwinger's effective Lagrangian EG] 



^Schwingcr — £~> \F ; Q) 



Here we have adopted the now common variables J18L h_ 
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The precise functional form of Schwinger's Lagrangian 
will not be needed: At the end of the calculation we shall 
see that it is sufficient to retain only the quartic terms 
beyond the Maxwell Lagrangian (the Euler-Heiscnberg 
Lagrangian pl|). 

Now in the Casimir geometry between parallel plates 
there is one additional invariant one could consider, 
namely 



H = n»F ua F° 



{E ■ nf - (fix Bf 



(4) 



where n M = (0, 0, 0, 1) is a unit vector orthogonal to the 
plates. (A similar invariant obtained by interchanging 
E and B is actually a redundant linear combination of 
the invariants T and TL.) In addition the coefficients in 
the effective Lagrangian can now explicitly depend on 
the z coordinate. When it comes to practical calcula- 
tions all these effects are suppressed by additional factors 



of H/(m e ca), but in the interests of generality we retain 
them for the time being and simply write: 



--effective 



C(FUx),x) 



(5) 



That is, the considerations of the following section are 
not limited to the Casimir parallel plate geometry but 
apply whenever the soft photon approximation (and the 
ancillary linearization procedure and restricted eikonal 
approximation) make sense. 



A. Equations of motion 

The complete equations of motion for nonlinear elec- 
trodynamics consist of the Bianchi identity, 



F[^,x] = 0, 
plus the dynamical equation 
/ dC 



= 0. 



(6) 



(7) 



We now adopt a linearization procedure: Split the elec- 
tromagnetic field into a background plus a propagating 
photon 



P _ Tp background _r_ f photon 

r 11V — r nil ~T J (U/ 



(8) 



Then assuming the background satisfies the equations of 
motion, and retaining only linear terms in the propagat- 
ing photon, we have 



p photon 
J[nu, A] 



0, 



and 




photon j 



0. 



dF, v dF a f 



(9) 



(10) 



(11) 



background 



equation ( JToj ) can be rewritten in the somewhat more 
compact form 



photon \ 



0. 



(12) 



Note that the tensor fl^ ual3 is symmetric with respect to 
exchange of the pairs of indices /iu and a/3, and antisym- 
metric with respect to exchange of indices within each 
pair. 

We now apply a restricted form of the eikonal approx- 
imation by introducing a slowly varying amplitude f^ v 
and a rapidly varying phase <p: 



£ photon 
J LLU 



f e i4 

J liv c 



(13) 
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The wave vector (actually a one-form, but we shall 
loosely refer to both vectors and one-forms as "vectors" 
in the text, because of the mapping provided by r]^ and 
its inverse rj^) is then defined as k^ — d^cf). This ap- 
proximation is similar to, but not quite identical with, 
the usual eikonal approximation. This is because one as- 
sumes that <j) varies on scales much smaller than those 
of the background, while, on the other hand, use of the 
Lagrangian (||) also implies that the components of k are 
much smaller than the values fixed by the electron mass 
(soft-photon regime). Under these hypotheses, 

Of™** k a f vP = 0. (14) 

But the background field is itself subject to quantum fluc- 
tuations, and to take this into account the coefficients of 
this equation are identified with the expectation value of 
the corresponding quantum operators in the background 
state \if}): 

^\n^\^)k a u P = Q. (15) 

In taking this expectation value we are using the fact that 
the fluctuations in the background fields are determined 
by the geometry (in the specific case of the Casimir geom- 
etry, by the distance between the plates), so that in the 
spirit of the restricted eikonal approximation there is a 
separation of scales between the background fluctuations 
and the propagating photon. 

The Bianchi identity (||) constrains f^ v to be of the 
form 

fnv = kf± a v — k v , (16) 

where we have introduced t he g auge potential a for the 
propagating field. Inserting (^6^ into (|l5|) we find 

{i/>\W m 'P\1>)k a kf l a v = 0. (17) 

In general, the tensorial quantity (ip\Q flat 'P\i{)) can be 
decomposed into an isotropic part plus anisotropic con- 
tributions, that we group together into a term A^ av>3 
with the same symmetries as f^""' 3 : 

(^Ifi^^l^) = rfi (if rf 3 - rf r! av ) + A^P, (18) 

where d\ is a function that can in principle be computed 
directly from the effective Lagrangian. 

B. Dispersion relations 

Equation ( |l7| ) represents a condition for a as a function 
of k — it constrains a to be an eigenvector, with zero 
eigenvalue, of the fc-dependent matrix 

A»»(k) = (V'lO"^^) k a k p . (19) 

Any non-zero solution corresponds to a physically pos- 
sible field polarization, that can be identified by a unit 



polarization vector e (provided a is not a null vector - 
a possibility that can always be avoided by a suitable 
gauge choice). 

A necessary and sufficient condition for the eigen- 
value problem A^ v a v = to have non-zero solutions 
is det (A^ v ) = 0; however, this gives us no information 
at all. Indeed, any a parallel to k is always a non-zero 
solution, so the condition det [A^ v ) = is actually an 
identity. On the other hand, a \\ k is merely an unphysi- 
cal gauge mode that corresponds to f^ u = by (fig), so 
we need to find other, physically meaningful, solutions of 
the eigenvalue problem. To this end, we exploit gauge 
invariance under a — > a + A k and fix a gauge, thus re- 
moving the spurious modes. It is particularly convenient 
to adopt the temporal gauge eto = 0. Then we can de- 

1/2 

fine a polarization vector e M = a^j (a v a v ) ' , and the 
eigenvalue problem A^ v e u — splits into the equation 

A 01 e t = 0, (20) 
plus the reduced eigenvalue problem 

A^e 3 =0. (21) 

The latter admits a nontrivial solution only if 

det (A ij ) = 0. (22) 

The condition (^) plays the same role as the Fresnel 
equation in crystal optics [E2| — it is a scalar equation 
for k and thus gives the dispersion relation for light prop- 
agating in our "medium" . 

An explicit calculation (see Appendix [^) gives 

det(A ij ) =oj 2 P 4 (k), (23) 

where to = —ka and V 4 (k) is a homogeneous fourth-order 
polynomial in the variables u> and fcj. This means that in 
the most general case there are four dispersion relations, 
corresponding to the four roots of the equation 

Vi(k) = 0. (24) 

But, if k = (—ui,k) is a root then so is — k = (ui,—k); 
thus, by CPT invariance, only two of these dispersion 
relations are physically distinct. 

Different polarization states are represented by linearly 
independent solutions of the eigenvalue problem (|2l]) , 
under the condition (|22|). [Obviously, ( po| ) cannot be 
independent of (|2l|), since we know that det(A^ v ) = 
0.] Thus, the space of polarizations is at most two- 
dimensional. Since ( [22] ) gives rise to two dispersion re- 
lations, the polarization states actually satisfy two (in 
general, different) eigenvalue equations, 

A^e^ = 0, (25) 

where r = 1, 2 labels the dispersion relations and A^ is 
obtained from A^ v by imposing the corresponding con- 
dition on k. Thus, in the general case, modes of the field 
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with different polarizations have different dispersion re- 
lations, hence propagate in different ways; this leads to 
the phenomenon of birefringence 

El Sill 



In some special cases (in particular, this behaviour is 
quite common in nonlinear electrodynamics) the polyno- 
mial Vi(k) factorizes into two quadratic forms, 



V 4 (k) = (7^ k^k v ) (7^ k a kp), 



(26) 



in which case we obtain two second-order dispersion re- 
lations: 



7(i) k/j,k u 



and 



T(2) kpk v 



0. 



(27) 



Effective metric 



We now want to consider the special situation in which 
not only does ( |26] ) hold, but also 7/^ = l^)- (That is, 
the fourth-order polynomial Vn(k) is a perfect square.) In 
this case one ends up with a single quadratic dispersion 
relation of the familiar form 



7 ^ k^ k v 



0, 



(28) 



where 7 MI/ is some symmetric tensor. It should be clear 
from our previous discussion that a necessary condition 
for this to happen is the absence of birefringence. Re- 
markably, we see that the wave vector is now null with re- 
spect to a (unique) "effective inverse metric" 7 ^ v . There- 
fore, the propagation of light can be described in terms 
of an effective geometry, defined by the metric tensor g^ v 
obtained by inverting 7 A " / , such that 7 A " / g J , p — S^ p . 

(Warning: Even when an effective metric is defined, we 
always raise and lower indices using the flat Minkowski 
metric r/^. In particular, note that g Mt/ ^ r\^ p r\ va ^ pa \ 
this justifies our use of different symbols, g and 7, for the 
effective metric and its inverse.) 

This situation implies that {ip\^l llo " J ^\'il)) must be alge- 
braically constructible solely in terms of 7^. In view 
of the symmetries of VL^ avl3 we know, without need for 
detailed calculation, that it must be of the form 



A (7^ 7 a/3 



(29) 



for some function A. 

Conversely, if VL^ au ^ is of the form (|2^) , then the ma- 
trix (|l!| is 

A^=A [-f^{^k a k )-(^ a k a )( 7 ^k )], (30) 

and the photon propagation equation A^ u a v = be- 
comes 

(7 aP k a k )j^ a, - (7 a a hp) J^K = 0. (31) 

This equation is obviously satisfied by the uninteresting 
gauge modes a \\ k, with no constraints on k. Solutions 
corresponding to a non- vanishing f UL , exist only if the co- 
efficient of / y fiv a v is zero, i.e., if (j28|) holds. Thus, the 



two polarization states propagate with the same disper- 
sion relation (|28), a nd there is no birefringence. 



Substituting (28) back into the propagation equation 
( |3l| ) wc find another relationship typical of this case, 



7"" k, L a v = 0. 



(32) 



Formally, the above equation looks like a gauge condi- 
tion. This might seem puzzling, because nowhere in the 
present subsection have we fixed a gauge. In fact, ( |3^ ) 
is a consequence of the dynamical equation A^ v a v = 0, 
when the "on-shell" condition ( |28| ) is satisfied, and it does 
not imply any gauge fixing. 

Finally, it is interesting to notice that there is now a 
self-consistency or "bootstrap" condition, 



3A7"" = (^|n^|^) ga0 . 



(33) 



This can be re-stated directly in terms of the fundamental 
coefficients as 



3 A 7"" 



di rf v (r) a0 gap) 
+A" ou " s g a . 



di rj' ia ga/3 V 



0u 



(34) 



We stress that these relations depend only on the as- 
sumed existence of a single unique effective metric g^ 
- they do not make any reference to other specifics of 
the quantum state. 



III. CASIMIR VACUUM 

Let us now consider a region of empty space delimited 
by two perfectly conducting parallel plates placed orthog- 
onal to the z axis at positions z = zq and z = zq + a, 
in the quantum state hp) — * \C) corresponding to the 
Casimir vacuum [Q, ||, g| . 



A. Symmetries and effective metric 

In the Casimir vacuum considerable information can be 
extracted by using only symmetry considerations, simi- 
larly to what Bryce DeWitt did for the stress-energy- 
momentum tensor in reference |2^|. Always working in 
the soft photon approximation, the presence of a pre- 
ferred direction and the symmetry of the configuration 
allow us to claim that the function d\ can only depend 
on the z coordinate, and to write A^^p in the form 

= d 2 {z) (5»» n a n - S"p n a n v 

+6 a p n v -6 a v n" n p ), (35) 

where d 2 {z) is another function. The calculation of d\(z) 
and d 2 (z) in this case is st raig htforward, using identities 
that follow from (|l8|) and (|3q): 



d 1 (z) = - (C \W\ U - 2 W\ v n M n a I C) ; 



(36) 
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d 2 (z) = - - {C\W\ V -4W av TV n a \ C) 



On denning the function 



and the tensor 



_ d 2 (z) 
~ di(z)' 



(C| r2 A£Q!l//3 1 C) takes the particularly simple form 

{C\W avt) \C) = d x {z) (7^ 7 Q/3 -7^ 7 va ). 



(37) 



(38) 



(39) 



(40) 



(pl|) . We find that e must be a solution of the problem 



Thus, from the discussion in subsection II C| , we deduce 
without further argument that there is no birefringence, 
and that 7^" given in ( |39| ) plays the role of an effective 
inverse metric for the propagation of light in the Casimir 
vacuum. The corresponding effective metric is 



(41) 



We stress that for this derivation of the effective metric 
to make sense it is necessary that the spatial variation in 
£(z) be compatible with the restricted eikonal approxi- 
mation: that is £(z) should vary slowly on the scale of 
the photon wavelength. This is certainly true for QED 
at lowest nontrivial order where we shall soon see that £ 
is in fact independent of z. 

The conclusion about the absence of birefringence can 
also be obtained more explicitly. Choosing the temporal 
gauge, and following the steps outlined in the general 
analysis of subsection II B, equation (EOT) becomes 



l^kiCj = 0, 

which also follows from (|32] ) when eo 
straightforward calculation gives 



(42) 



0. A tedious but 



det(A ij ) = -(1 + df 



(7^ K 



(43) 



Thus the dispersion relation determined from (|22|) takes 
in this case the form (^8|), as expected. 

We note here that, although an equation of the form 
( p8| ) has been already obtained by Dittrich and Gies jl8) , 
the expression for the effective metric underlying their 
dispersion relation is subtly at variance with our results 
( |39| ) and (JO]). We shall comment about the reason for 
this discrepancy in Appendix 



B. Polarization states 

In order to explicitly find the polarization states, let us 
evaluate jl9| ) using the expression for (C\£W' av P\C) and 
(Eq), and then consider the reduced eigenvalue problem 



T 3 u 



0, where 



.4' 



(44) 



This eigenvalue problem, however, is manifestly equiv- 
alent to the single equation ([l2|), which is thus the 
only constraint that the polarization states must sat- 
isfy. Therefore, the space of polarization states is two- 
dimensional, as expected. 

A basis for such a space can be easily constructed by 
considering generic propagation in the xz plane [not a 
restrictive hypothesis, because of 0{2) invariance with 
respect to rotations around the z axis] , described by the 
3-vector (these are taken to be covariant components, 
index down): 



k = (jfc|sin6>,0, fc| cos 6 1 

Equation ( fi"2| ) then becomes 

sin 6 ei + (1 + f ) cos 6 e 3 = 0, 
so two independent polarizations are: 
e« = (0,1,0); 



(45) 



(46) 



'(2) _ 



N{6) 



((l + £)cos0,O,-sin0), (47) 



where 



N{0) = (1 + 2 £ cos 2 + £ 2 cos 2 6) 1/2 



(48) 



is a normalization coefficient. Note in particular that 
e *- 2 ) is not perpendicular to k when viewed in terms of 
the Minkowski metric rf LV , though they are perpendic- 
ular when viewed in terms of the effective metric 7^". 
Furthermore we have chosen to make a unit vector 
with respect to the Minkowski metric, not with respect 
to the effective metric. 



C. Phase, signal, and group velocities 



If k has the form (|45|), equation (|2£]) becomes (the 
norms are with respect to the Euclidean spatial metric 
induced by rj^): 



oj 2 = \k\ 2 +i\k\ 2 cos 2 ( 
The phase velocity is given by 



^phasc(0) 



to 

W\ 



(1 + £ cos 2 6 



1/2 



(49) 



(50) 



and is independent of the polarization. This is again a 
consequence of the fact that the quantity 7 ^ fc M k v ap- 
pears squared in det (j4 y ), so equation ( p2| ) describes a 
degenerate fourth-degree surface in the space of the vec- 
tors k/ui. Hence, there is only one dispersion relation, 
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equation (|49|), independent of the polarization state, and 
only one phase speed for each value of |fc|. This confirms 
that birefringence does not take place in the Casimir vac- 
uum. 

Note that if £ > and ^ tt/2, we have v p hase(0) > 1. 
Since, in the limit co — > +oo, v p hase equals the front veloc- 
ity (signal velocity) j7j , one is tempted to argue that the 
propagation is supcrluminal for all values of different 
from 7r/2. This conclusion can also be inferred directly 
from ( p8[ ) , which for £ > implies that k is timelike (with 
respect to the undisturbed Minkowski metric). Since k 
can also be interpreted as the four- vector orthogonal to a 
surface of discontinuity of the field Jll| [l7) , it follows that 
such a surface is spacelike (with respect to the Minkowski 
metric), i.e., that electromagnetic signals travel "faster 
than light" . (This phrasing is standard but unfortunate, 
and is logically indefensible. It should always be inter- 
preted in the sense "faster than light would have travelled 
in an undisturbed portion of normal vacuum".) Equiva- 
lently, ( fll"| ) implies that for £ > the lightcones of the ef- 
fective metric g^ v are wider than those of the Minkowski 
metric f]^- Unfortunately, while it is certainly true that 
our treatment is valid at high frequencies (with respect to 
those associated to the background scales), it neverthe- 
less also requires u) <C m c — the condition under which 
one can use the Lagrangian ([!]) — , so we have no direct 
information about the strict to — > +oo limit. See, how- 
ever, reference 0] for an indirect argument based on the 
Kramers-Kronig dispersion relation that combined with 
the present calculation is sufficient to establish superlu- 
minal propagation for the signal velocity. In the special 
case 9 — tt/2 one has rj^k^k^ — 0: Photons propagat- 
ing parallel to the plates travel at the standard speed of 
light. 

The group velocity is a little tricky, it equals 



dui 



1 



£(fc 



(51) 



so equations ([15]) and ( pOj ) give 
1 



^group (^) 



^phasc (0) 



(sin0,O,(l+£)cos0). (52) 



In particular, the group velocity v gTonp is not parallel to 
the wave-vector k, though it is always orthogonal to the 
polarization vector. Taking the norm (in the physical 
Minkowski metric), 



^group 



N(9) 



1 + (2£ + £ 2 )cos 2 



1 + £ COS 2 



1/2 



(53) 



The group speed is at all angles slightly greater than (or 
at worst equal to) the phase speed. (See figure 1.) Indeed 



group 



(0) 



V 



phase 



(0) 



£ 2 cos 2 sin 2 



phase 



(0) 



(54) 



(See figure 2.) Care should be taken to realize that here 
is the angle between the wave vector and the normal, 



and is not quite the same as the direction of propaga- 
tion of the wave packet. Fortunately at both normal 
incidence (0 = 0) and parallel propagation (0 = tt/2) 
the distinction between group and phase velocities dis- 
appears. Furthermore, as we shall see below, if we work 
to second order in the fine structure constant, the differ- 
ence between group and phase speed is negligible at all 
angles, although a difference of order a 2 still remains be- 



tween the directions of - 



'group 



and k. The group velocity 



( p3| ) is greater than 1 when £ > 0. Thus, the group speed 
is larger than 1 whenever the phase speed is, although in 
general their individual values are different. 



1.2 1.4 



FIG. 1: Group and phase speeds as a function of angle, 
from 6 — to 9 — n /2 radians. (The group speed is the 
upper curve.) For clarity we have greatly exaggerated the 
physically expected value of £ by setting £ = 1/2. 



IV. SIZE OF THE EFFECT 

So far we have seen that, based solely on symmetry 
considerations, we can eliminate the possibility of bire- 
fringence in the Casimir vacuum and place strong con- 
straints on the general features of photon propagation 
in this vacuum. Our results are in fact generic to any 
form of nonlinear electrodynamics subject to the bound- 
ary conditions appropriate for the Casimir vacuum, and 
are not specifically restricted to QED. Where QED is 
important is in determining the specific form of the func- 
tions d\(z) and d,2(z), which determine £(z), and hence 
determine the size (but not the qualitative features) of 
the Scharnhorst effect. 

Unless one wants to perform calculations to orders 
higher than a 2 , it is sufficient to consider the Euler- 
Heisenberg Lagrangian plfl which, in the T-Q formalism 
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C 2 G 




FIG. 2: Difference between the squares of the group and 
phase speeds as a function of angle, from 6 = to 9 = 7r/2 ra- 
dians. For clarity we have greatly exaggerated the physically 
expected value of £ by setting £ = 1/2. 



adopted above, takes the form 



C 



EH 



1 

47T 



jF+c a T 2 + c 2 g 2 , 



with 



Cl 



cr 



907r 2 m 4 ' 



C2 



7a 2 



3607T 2 TO 4 



(55) 



(56) 



The terms proportional to T 2 and Q 2 of this Lagrangian 
are quartic in the field, and describe the low-energy limit 
of the box diagram in QED, when four photons couple 
to a single virtual electron loop. Thus, the Lagrangian 
(|55|) is only accurate to order a 2 , and it is meaningless 
to retain higher order terms within this model. 

In particular, deviations from (3+1) Lorentz symme- 
try due to the presence of the plates (the plates re- 
duce the symmetry group to that of (2+1) Lorentz in- 
variance) must on physical grounds vanish as the plate 
separation goes to infinity, where one must recover the 
full (3+1) Lorentz symmetry. On dimensional grounds 
this implies that such terms will be suppressed by some 
function of h/ (m e ca) and can, to lowest nontrivial order, 
be neglected in the effective Lagrangian. On the other 
hand, even for lowest nontrivial order in the effective La- 
grangian (that is, for the Euler-Heisenberg Lagrangian), 
the presence of the plates leads to a nontrivial expecta- 
tion value for (ClF^ F a p\C) and will in this way con- 
tribute to the effective metric. 

For the Euler-Heisenberg Lagrangian, the tensor 
W a p is 



fin 

E a/3 



+ £L f" v F a/3 + ^ *F» V *F t 



a/3, 



so in the Casimir vacuum 



dl(z) = "l6. +0(a2) 



(57) 



(58) 



while d 2 (z) is of order a 2 . Hence, to first order in a 2 , we 
have £(2) = — I6ird 2 (z), and 



-Uphasc(#) = Ugroup(0) = 1 - 8nd 2 (z) COS 2 9. 



(59) 



Though in principle the coefficient d%{z) could depend on 
z, the position relative to the two plates, we shall see that 
in the specific case of the Casimir vacuum it is simply a 
position-independent number. 

To establish this, start with (|37|), and insert the specific 
form (HI) ■ Using the algebraic identity Jl|, ||(| 



* ppv * p 
one easily gets 



F F av 2 J~ S^ a , 



(60) 



+2tt ( Cl + c 2 ) T\ + ( Cl ~ C2 )- F S \, (61) 



where T^ a is Maxwell's stress-energy-momentum tensor. 
Performing the indicated traces, and using the fact that 
the Maxwell stress-energy tensor is traceless, we find 

d 2 (z) = ^-( Cl +c 2 ) (C\T ZZ \C). (62) 

The symmetries of the Casimir vacuum stress-energy (as 
analyzed for instance by DeWitt then imply 



d 2 (z) = 4tt(ci + ca) (C\T 00 \C). 

Finally the well-known result (C\T q\C) = 
allows us to write 



d 2 = 



llTra 2 



64800 a 4 mi 



(63) 
-7r 2 /720a 4 

(64) 



And in particular, this implies d 2 (z) is position indepen- 
dent. Thus, at first order in a 2 , 

ll7T 2 a; 2 

VphascW = Voup(0) = 1 + 8 i 00q 4 m 4 C0S2 6 - ( 65 ) 

This expression reproduces Scharnhorst's result in the 
case = 0, generalizing it to an arbitrary direction of 
propagation. It is interesting to notice that the correction 
is essentially determined by the expectation value of the 
energy density, (C\Tqq\C), in agreement with the general 
results of IO, Ell. 
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V. CONCLUSIONS 

We have outlined a general scheme that allows one 
to write down the dispersion relation, and to find the 
polarization states, for electromagnetic radiation propa- 
gating in a region where nonlinear effects cannot be ig- 
nored. The deviations from the behaviour in Maxwell's 
theory are completely described by the tensorial quan- 
tity A''" 1 '' 3 , that has its origin in the anisotropy of the 
medium. In particular, we have investigated the case 
of propagation in the Casimir vacuum, where there is a 
privileged direction in space, identified by a unit vector 
n. We have seen that symmetry considerations alone im- 
ply that lS. p - a v must be of the form (|35|). This, in turn, 
implies that there is only one value for the speed of light, 
independent of the polarization. Thus, the possibility of 
birefringence in the Casimir vacuum is completely ruled 
out by very general arguments (the abstract form of the 
Lagrangian for nonlinear electrodynamics, plus the ten- 
sor structure of t^ a v p dictated by the geometry). Be- 
cause of its generality, this conclusion applies not just to 
QED itself but to arbitrary types of nonlinear electrody- 
namics (such as for instance, Born-Infcld theories). 

How are we to interpret this result in view of the 
fact that nonlinear electrodynamics generically does lead 
to birefringence Jl5|, [H], [lTj? The key point is that in 
those analyses the background field (-F^background is 
nonzero. More generally, even after averaging over quan- 
tum fluctuations of the background, (iplF^lip) 7^ in 
the quantum state appropriate to those analyses. In the 
Casimir vacuum on the other hand, the expectation val- 
ues linear in the field do vanish, {C\F^ V \C) = 0, and 
it is only the quadratic expectation values (and higher) 
that are nonzero ({C\W af} \C) ^ 0Y This key differ- 
ence makes the analyses of jLq, |l^, [17| inapplicable to 
the Casimir vacuum — technically speaking, the polariza- 
tion basis used in those papers fails to be meaningful 
for the Casimir vacuum. On the other hand, the discus- 
sion of the present paper could be easily adapted to treat 
light propagation in the presence of an electromagnetic 
background field. It is sufficient to keep in mind that in 
this case there is a preferred 2-dimensional plane defined 
by the 2-form (i^lF^lifj) . Correspondingly, the form of 
w jvj kg more complicated [see (|A2|)]. 

We reiterate that the absence of birefringence is crucial 
for the use of the "effective geometry" approach (adopted 
in this type of context by Latorre et al p"2| ], and fur- 
ther developed by Novello et al |L6| , and by Dittrich and 
Gies |l£| ) . Only if the propagation of light does not de- 
pend on its polarization and is thus, in a sense, univer- 
sal, it is meaningful to describe it by a single effective 
metric. As a consequence of our analysis, photon prop- 
agation in the Casimir vacuum can indeed be phrased 
entirely in terms of the effective metric (39). This obser- 
vation is potentially important in that "effective metric" 
approaches similar in spirit to the above are currently at- 
tracting attention in fields as diverse as acoustics |24|, , 
optics p(| p7), superfluid quasiparticles 28 , and Bose- 



Einstein condensates p9[ . 

Finally, we mention related work on the light cone con- 
dition for a thermalized QED vacuum due to Gies |3^| 
(wherein the analysis implicitly relies on taking the quan- 
tum expectation value of VL^ ua " in a manner somewhat 
analogous to the present paper), and intriguing results 
on photon propagation in rather general linear theo- 
ries in classical backgrounds due to Obukhov and co- 
workers fn], I32I, I33I1 . 
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APPENDIX A: CONNECTION WITH OTHER 
FORMULATIONS 

We have tried to set up the formalism in a streamlined 
and self-contained manner. Nevertheless to aid compar- 
ison with other results in the literature it is useful to 
give the explicit form of the tensor for generic 

Schwinger-like Lagrangians of the form (|l|) . For such La- 
grangians 



-- (dg£.) t^ aP 



M^ Q/5 , (Al) 



where 

M »»af3 = ppv pa/3 (g2 + * p^ * pa(3 (q2j^ 

+ (pi* *p<*P + *p^ pot?} d:/rgC ( A2 ) 

has the same symmetries as W a P. 

As soon as one inserts this tensor into the photon equa- 
tion of motion (14), the completely antisymmetric part 
proportional to the Levi-Civita tensor drops out, because 
of the Bianchi identity @. The remaining pieces repro- 
duce the photon equation of motion in the perhaps more 
usual form considered by Dittrich and Gies jl8), or Nov- 
ello and co-workers [|l6| p"7| |. 

Also note that, depending on the details of the geom- 
etry and the quantum state, M^ va $ can contribute to 
both di and to A^ Il/Q ' 3 . 



APPENDIX B: FRESNEL EQUATION IN 
NONLINEAR ELECTRODYNAMICS 



Our purpose is to prove equation (23). The spatial 
components of the matrix (|l9|) are 

(Bl) 



+k m k n (v>|n* m ' n |v> 
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Let us define the unit vector k = k/\k\. Then the com- 
ponents of A l i in a basis with one axis directed along k 
are 

A ij kj = u? (ip\n i0j0 \ip) kj-uj k m (ip\{t mj0 \ip) kj = wv\ 

(B2) 

In particular 

A^kikj = to 2 (tp\n iOj0 \tp) kk = lu 2 S. (B3) 
Then the matrix A lJ has the following structure: 



LU 2 S L0V J 

toV 1 T IJ 



(B4) 



where / and J label the two directions orthogonal to k. 
Evaluating the determinant by expanding in the first row 
or column, it is easy to see that every term will contain 
at least two factors of w, which establishes equation ( |23| ) 
as desired. 



Specializing to the Casimir vacuum, and normalizing ap- 
propriately one has 



2d 2 



2d x 



(C6) 



clearly different from (p9|). Nevertheless, when one con- 
siders QED corrections to the Maxwell Lagrangian, d\ is 
of order a , while d 2 is of order a 2 (see section above); 
therefore, the two metrics (|9|) and (|C^ ) agree to order 
a 2 . 

The reas on f or the discrepancy is that the formal re- 
placement ( |C3| ) is limited in a subtle manner. Strictly 
speaking, it is correct as it stands only if photon propa- 
gation is governed by the Minkowski metric rj^; which 
is exactly the situation we are trying to get away from. 
In order to see this explicitly, let us introduce an or- 
thonormal tetrad (t^em, ep)) w ) > with u timelike and w 
spacelike, so one can write 



APPENDIX C: POLARIZATION SUM 



E 



(C7) 



In this appendix we discuss the reason for the difference 
between our expressions ( |39| ) and (|4l]) for the effective 
metric and those implicit in the "light cone condition" 
in reference (il|. In the spirit of that paper a dispersion 
relation would be derived as follows (a certain amount of 
"translation" is required to get that formalism to match 
with the current notation). First, one writes equation 
( |l7| ) for two polarization states and orthogonal 
to each other: 

{^|fi^|V) k a k e^ = o. (ci) 

Working in Lorentz gauge k ■ = this can be written 
as 

d x k 2 e/t, + A» a „ k a k^y, = 0, (C2) 



where e/\ = rf v t^. One then takes the scalar prod- 
uct of each equation with the corresponding polarization 
vector. Finally, the resulting two equations are summed. 
Furthermore, it was explicitly asserted that one could 
effectively replace |13] 



E^W 



(C3) 



in the terms containing A^^p. Under this hypothesis 
one obtains, following the steps described above, 



2d l k 2 + A» a ua K k v =0, 



(C4) 



which is equivalent to equation (14) of reference [g_8|. 
This has the form of the dispersion relation 7 ^ k^ k v — 
0, with 



7"" oc 



2 di if + A^^ n a 



Va/3- 



(C5) 



Working in the Lorentz gauge, as in reference |18 , the 
two-dimensional subspace of Minkowski spacetime, or- 
thogonal to ?(!) and ?(2) and spanned by the unit vectors 
u and 10, contains the wave vector k, so 



(k -u) it M + (k ■ w) w M . 



(C8) 



Following the same steps that led to (|C4j), but using (C7) 
instead of the replacement (C3), we find, in place of the 
single term A^^ k^ k v occurring in (C4), the two terms 



A^ a va k^k v + k 2 A liavP u^u"w a w f) 1 (C9) 



where we have used (C8) and the symmetry properti es o f 
Atuxup- Clearly, the second term in the expression (C9) 
cannot be ignored unless k 2 — 0. We conclude that it is 
the subtly incorre ct re pl acem ent ( p3|) which implies that 
the derivation of (|C4]), (|C5|), and (|C6|) is limited to first 
order in a 2 . For the Casimir vacuum, 

A^avfj u v w a w f3 = d 2 [(n ■ u) 2 - {n ■ w) 2 ] 



2 

E 

r=l 



(n ■ < 



■ Hr)) 



(CIO) 



Using now the relation k 2 (n-Zf r )) = 0, which follows 
from ( |C2| ) upon contraction with n in the Lorentz gauge, 
we see that 



k A lJ , av/ 3U fl u u w a w 



-da k 2 , 



(Cll) 



so (C4), corrected by the additional term ( C10 ), indeed 
reproduces the metric (|39|). 

In general, if the photon dispersion is in fact deter- 
mined by some (unique) effective metric g^ v then one 
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can choose a tetrad like fit, em 
thonormal with respect to g^,,, 



so now 



but which is or- 



E 

r=l 



£ (r) £ (r) 



7"" + ^ U" 



(C12) 



Applying this polarization sum to the photon propaga- 
tion equation written in the form (|l7j), and using the 
properties ( p8| ) and (|3^), we see that the strictly correct 
replacement for the polarization sum is 



E ? A 



(r) 



(C13) 



This modified polarization sum, applied to (|l7|), is con- 
sistent with the "bootstrap" condition (|33[). 

In the particular case of the Casimir vacuum, it is a 
simple matter to verify, using our expression for 7 that 
the bootstrap condition is indeed satisfied. It is now also 
clear why the formalism of [jl8| works to first order in a 2 : 
Since the difference between rf 1 ' and 7 M1/ is itself of order 
£ (that is, of order a 2 ) the difference between t he (s ubtly 



incorrect) formula (C6), derived on the basis of (C3), and 
the consistency condition (^) is automatically of order 
£ 2 (a 4 ). This may also be explicitly verified using the 
approximation £ <C 1 so that 



(C14) 



Substituting this into ( p3[ ) self-consistently reproduces 
7 M!y to the required order. Note that in almost all cases 
of interest one is quite content to work to first order in 
a 2 in which case all of these subtleties are moot. 

We conclude with some general comments about the 
derivation of a metric from a polarization sum: If there 
is only one dispersion relation for all polarizations, as for 
instance in the Casimir vacuum discussed above, then 
e^ 1 -* and ?^ 2 ^ both satisfy the same equation 



A^(k) ?W = 



(C15) 



[see (p5|)], with the same matrix A^ v (k). However, if 
e^ 1 -* and correspond to two different dispersion re- 
lations, their propagation equations will actually differ, 
because now ^ Ar 2 \ . Therefore, a procedure of the 
form described at the beginning of this appendix is fun- 
damentally flawed in the general case, as it fails to take 
into account the possibility of birefringence, and is in 
fact (strictly speaking) incompatible with birefringence. 
Indeed, summing the equations 



(i) 



and 



-TV.U (2) (2) _ 
A (2) e n e v — U ; 



(C16) 



which do take birefringence into account, appears rather 
impractical, and it is not obvious whether the structure 
of the resulting equation will allow the use of a polar- 
ization sum at all. Thus, producing "average dispersion 
relations" by polarization sums is at some deep level in- 
ternally inconsistent, except in the cases when no polar- 
ization sum is needed. It is only if one is limiting interest 
to lowest-order (a 2 ) corrections away from Minkowski 
space that the polarization sum makes sense for a bire- 
fringent system, and only then in the sense that 



E 



€ (r) £ (r) 



0(a 2 ). 



(C17) 



That is: Although equation (|17| ) is exact, we could choose 
to contract it with approximate polarization states appro- 
priate to flat Minkowski space. As long as the deviations 
from the ordinary speed of light are small, so are the devi- 
ations of the true polarization states from the usual ones, 
and then an approximate polarization sum of the above 
form is useful. We conclude that the polarization sum 
technique is most useful if there is no birefringence, and 
that in the presence of birefringence it is, at best, only 
useful when strictly limited to lowest-order corrections to 
Minkowski space propagation. 



[1] 

[2] 

[3] 
[4] 
[5] 

[6] 



H. B. Casimir, "On the attraction between two per- 
fectly conducting plates," Kon. Ned. Akad. Weten- 
sch. Proc. 51. 793-795 (1948). 



V. M. Mostepanenko and N. N. Trunov, The Casimir 



[7] 



Effect and its Applications (Oxford Science Publications, [8] 
Clarendon Press, Oxford, 1997). 

G. Plunien, B. Muller and W. Greiner, "The Casimir 
effect," Phys. Rep. 134, 87-193 (1986). 

K. Scharnhorst, "On propagation of light in the vacuum [9] 

between plates," Phys. Lett. B236, 354-359 (1990). 

G. Barton, "Faster-than-c light between parallel mirrors. 

The Scharnhorst effect rederived," Phys. Lett. B237, [10] 

559-562 (1990). 

G. Barton and K. Scharnhorst, "QED between parallel 
mirrors: light signals faster than c, or amplified by the [11] 



vacuum," J. Phys. A 26, 2037-2046 (1993) 
K. Scharnhorst, "The velocities of light in modified QED 
vacua." An n. Phys. (Leipzig) 7, 700-709 (1998) ]hep ' 
th/981022l| 



I. T. Drummond and S. J. Hathrell, "QED vacuum polar- 
ization in a background gravitational field and its effect 
on the velocity of photons," Phys. Rev. D 22, 343-355 
(1980). 

R. D. Daniels and G. M. Shore, " 'Faster than light' 
photons and ch arged black hole s," Nucl. Phys. B425, 
634-650 (1994) ]hep-th/9310i"l4| . 

R. D. Daniels and G. M. Shore, " 'Faster than light' 
photons and rotating black holes," Phys. Lett. B367, 
75-83 (1996) fer-qc/950804{| . 

G. M. Shore, " 'Faster than light' photons in gravitational 



11 



[12] 

[13] 
[14] 
[15] 



[16] 
[17] 

[18] 

[19] 
[20] 
[21] 
[22] 



[24] 



fields — Causality, ano malies and hori zons." Nucl. Phys. 

B460, 379-394 (1996) |gr-qc/950404l| . 

J. I. Latorre, P. Pascual and R. Tarrach, "Speed of light 

in non-trivial vacua," Nucl. Phys. B437, 60-82 (1995) 

flrep-th/9408016 |. 

A. D. Jackson, A. Lan de and B. Lautm p, "Apparent 
superluminal behavior," physics/0009055] . 
A. M. Steinberg, "No thing goes faster than light," Phys. 
World 13 (9), 21 (2000). 

S. L. Adler, "Photon splitting and photon dispersion in 
a strong magnetic field," Ann. Phys. (N.Y.) 67, 599-647 
(1971); 

"Comment on "Ph oton splitting in st rongly magnetized 
objects revisited" ," ta,stro-ph/9601156i 
S. L. Adler and C. Schubert, "Photon splitting in a strong 
magnetic field: Recalculation and comparison with previ- 
ous calculations," Phys. Rev. Lett. 77, 1695-1698 (1996) 



(1998) |gr-qc/97120ic| ]: 

"Hawking radiation without black hole entropy," 
Phys. Rev. Let t. 80, 3436-3439 (1998) 



gr-qc/9712016 



[26] 



"Acoustic black holes," |gr-qc/9901047 

S. Liberati, S. Sonego and M. Visser, "Unexpectedly large 

surface gravities for acoustic horizons?" Clas s. Quantum 

Grav. 17, 2903-2923 (2000) ]gr-qc/0003105| . 

U. Leonhardt and P. Piwnicki, "Optics of nonuniformly 

moving media," Phys. Rev. A 60, 4301-4312 (1999) 

[ physics/9906038} ; 

"Relativistic effects of light in moving media with ex- 
tremely tow_grou_D_ve]acity," Phys. Rev. Lett. 84, 822-825 



[27] 



cond-ma 



t/9906332| . 



| hep-th/9605035 | 
M. Novello, V. A. De Lorenci, J. M. Salim and R. Klip- 
pert, "Geometrical aspects of light propagation in non- 
linear electrody namics," Phys. Rev. D 61, 045001 (2000) 
fer-qc/9911085 |. 

V. A" De Lorenci, R. Klippert, M. Novello and 
J. M. Salim, "Light propagation in non-linear el ec- 
trodynamic s," Phys. Lett. B482, 134-140 (2000) 



(2000) _ 

M. Visser, "Comment on Relativistic effects of light 
in moving media with extremely low group ve locity" , 
Phys. Rev. Lett. 85, 5252 (2000) fer-qc/000201~ll ]. 
U. Leonhardt and P. Piwnicki, "Reply to the Comment 
on Relativistic Effects of Light in Moving Media with Ex- 
tremely Low_Gwu_pVelocity" , Phys. Rev. Lett. 85, 5253 



[28] 



(2000) 

G. E. Volovik 



gr- 



qc/00030ie| ]. 



"Simulation of quantum field theory and 



gravity in superfluid He," Low Temp. Phys. (Kharkov) 



qc/0005049| 



24, 127-129 (1998) ]cond-mat/9706172j ] 



ial QED vacua," Phys. Rev. D 58, 025004 (1998) I 


aep- 


Zh. Eksp. r 


ph/9804375]. 




mat/9712187 



N. B. Kopnin and G. E. 
locity and event horizon in 
with "relativistic" fermionic 
Fiz. 67, 



W. Dittrich and H. Gies, Probing the Quantum Vacuum, 
Springer Tracts in Modern Physics 166 (Springer, 2000). 
J. Schwinger, "On gauge invariance and vacuum polar- 
ization," Phys. Rev. 82, 664-679 (1951). 
W. Heisenberg and H. Euler, "Consequences of Dirac's 
theory of positrons," Z. Phys. 98, 714-732 (1936). 
L. D. Landau, E. M. Lifshitz and L. P. Pitaevskii, Elec- 
trodynamics of Continuous Media, 2nd edition (Oxford, 
Pergamon Press, 1984). 



Volovik, "Critical ve- 
pair-correlated systems 
quasiparticles," Pisma 
124-129 (1998) flcond- 



G. E. Volovik, "Gravity of monopole and string and 
gravitational constant in 3 He-A," Pisma Zh. Eksp. Teor. 
Fiz. 67, 666-671 (1 998); JETP Lett. 67, 698-704 (1998) 
|x3nd-mat/980407g ]. 



[29] 



G. E. Volovik, "Li nks between gr avity and dynamics of 
quantum liquids," ^r-qc/0004049 

L. J. Garay, J. R. Anglin, J. I. Cirac and P. Zoller, "Sonic 
analog of gravitational black holes in Bose- Einstein con 



[23] B. S. DeWitt, "Quantum gravity: the new synthesis," in 



densates," P hys. Rev. Lett. 85, 4643-4647 (2000) [?r- 



qc/0002015 



General Relativity, edited by S. W. Hawking and W. Is- 
rael (Cambridge, Cambridge University Press, 1979), 
pp. 680-745. 

W. G. Unruh, "Experimental black hole evaporation?" 



"Sonic black holes in dilute Bose -Einstein conde nsates 
Phys. Rev. A 63, 023611 (2001) |gr-qc/0005131 |. 
[30] H. Gies, "Light cone condition for a thermalized 
QED vacuu m," Phys. Rev. D 60, 105033 (1999) |iep 



Phys. Rev. Lett. 46, 1351-1353 (1981); 



"Dumb holes and the effects of high frequencies on Mack [3lT 
hole evaporatio n," Phys. Rev. D 51, 2827-2838 (1995), 
fer-qc/9409008 |. 

(Title changed in journal: "Sonic analog of black holes [32] 
and...") 



ph/9906303] 



Y. N. Obukhov and F. W. 
from li near electrodyna mics," 



[25] M. Visser, "Acoustic propagation in fluids: An unex- 



pcctcd exampl e of Lorentzian geometry," 
gr-qc/9311028| ; 

"Acoustic black holes: Horizons, ergospheres, and Hawk- 
ing radiation," Class. Quantum Grav. 15, 1767-1791 



Hehl, "Spacetime metric 
Phys. Lett. B458, 466-470 
(1999) |gr-qc/9904067 |. 

F. W. Hehl, Y. N. Obukhov and G. F. Rubilar, 
"Spacetime metric from linear electrodynamics. II, 
qc/9911096f 

Y. N. Obukhov, T. Fukui and G. F. Rubilar, "Wave 
propagation in linear electrody namics," Phys. Rev. D 62, 
044050 (2000) |gr-qc/0005018|. 



